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ABSTRACT 

A simple function for the spectral power P(e,i) = vL(v) is proposed to model, with 9 
qq ' parameters, the spectral and temporal evolution of the observed nonthermal synchrotron power 

■ flux from GRBs in the blast wave model. Here e = hv/m e c 2 is the observed dimensionless 
photon energy and t is the observing time. Assumptions and an issue of lack of self-consistency 
are spelled out. The spectra are found to be most sensitive to the baryon loading, expressed 
in terms of the initial bulk Lorentz factor To, and an equipartition term q which is assumed to 

Xf^ ' be constant in time and independent of Tq. Expressions are given for the peak spectral power 

P p (t) — P(e p ,t) at the photon energy e = e p (t) of the spectral power peak. A general rule is 
that the total fireball particle kinetic energy Eq ~ Hotd, where td oc T 8 ^ 3 is the deceleration 
time scale and IIo = P(e p , td) oc T^ 3 is the maximum measured bolometric power output in 
' radiation, during which it is carried primarily by photons with energy £q — e p (td) oc qTg. 

t~^> . This rule governs the general behavior of fireballs with different baryon loading. Clean 

fireballs with small baryon loading (ro 3> 300) are intense, subsecond, medium-to-high energy 

■ 7-ray events, and are difficult to detect because of deadtime and sensitivity limitations of 
previous 7-ray detectors such as EGRET on CGRO. Dirty fireballs with large baryon loading 
(ro <C 300) produce transient emissions which are longer lasting and most luminous at X-ray 
energies and below, but these events are lost behind the glow of the X-ray and lower-energy 
background radiations except for rare serendipitous detections by pointed instruments. The 
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correlation between hardness and duration of loaded GRB fireballs (100 < T < fOOO) follows 
,1— > . from this rule. 
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1. Introduction 



5^ 

The cosmological origin of GRBs has been established as a result of follow-up observations of fading 
X-ray counterparts to GRB sources discovered with the Beppo-SAX mission (e.g., Costa et al. 1997; Metzger 
et al. 1997; van Paradijs et al. 1997; Djorgovski et al. 1997; Kulkarni et al. 1998a). The time profiles of 
the X-ray afterglow light curves are generally well fit by power laws with temporal indices x m the range 
1.1 < X < 1.5 where, at a fixed observing frequency, the flux tfi oc t~ x (e.g, Costa et al. 1997; Feroci et al. 
1998). Only the X-ray afterglow behavior of GRB 970508 is not well fit with a single power law, as its X-ray 
flux displays a peak between 1-2 days after the GRB. Its overall X-ray behavior from 30 seconds to 1 week 
following the burst event does, however, follow a power-law decay with ^ = 1.1 (Piro et al. 1998a). The 
light curves of GRB counterparts detected at optical and radio frequencies also generally display power-law 
behavior with a similar range of temporal indices (e.g., Sahu et al. 1997; Galama et al. 1998; but note Groot 
et al. 1998; Frail 1998). 
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The long wavelength afterglows of GRBs are most simply explained with the fireball/blast-wave model, 
and were predicted prior to their discovery (Paczyfiski & Rhoads 1993; Meszaros & Rees 1997). In this 
model, the impulsive release of a large amount of energy in a small volume is transformed into relativistic 
plasma whose baryon-loading parameter, expressed in terms of the initial bulk Lorentz factor T of the 
fireball after it has become optically thin, obtains values much greater than unity (Blandford & McKee 
1976; Cavallo & Rees 1978; Shcmi & Piran 1990; Rees & Meszaros 1992; Piran & Shcmi 1993). The 
dominant radiation mechanism producing the prompt GRB emission and radio through X-ray afterglows is 
nonthermal synchrotron emission (Meszaros & Rees 1993a; Katz 1994; Tavani 1996; Cohen ct al. 1997). The 
power-law behavior of the long wavelength afterglows is a consequence of the deceleration of the blast wave 
as it sweeps up material from the ambient medium (Wijers, Rees, & Meszaros 1997; Waxman 1997a; Vietri 
1997a; Dermer & Chiang 1998). The blast-wave model also provides an explanation for the > 100 MeV 
emission observed with the EGRET instrument from bright BATSE bursts (Dingus 1995; Catelli, Dingus, 
& Schneid 1996; Hurley et al. 1994) as due cither to synchrotron self-Compton (SSC) emission (Meszaros, 
Rees, & Papathanassiou 1994; Chiang & Dermer 1998) or to the nonthermal proton synchrotron radiation 
from hadrons accelerated in the relativistic blast wave (Vietri 1997b; Bottcher & Dermer 1998). 

In this paper, we propose a simple function to describe the observed prompt and delayed emissions 
from a decelerating blast wave in the fireball/blast- wave model. We use this function to perform a 
parameter study, and find that the temporal and spectral behavior of the blast-wave emissions are robust 
to orders-of-magnitude changes in all parameter values except for the baryon-loading parameter Tq and an 
cquipartition parameter q. If q is not strongly dependent on r , as quantified below, then the observability 
of a fireball is most strongly determined by the value of the baryon-loading factor To. Clean fireballs with 
small baryon loading (r » 300) produce intense, subsecond, mcdium-to-high energy (» 10 MeV) 7-ray 
events, and heavily-loaded dirty fireballs (ro <C 300) produce longer duration transients which are most 
luminous at X-ray energies and below. Limitations of earlier telescopes may have prevented the discovery of 
clean fireballs, though some previous serendipitous detections of X-ray transients might represent members 
of the dirty fireball class. 

In Section 2, we introduce a function to model the temporally-evolving nonthermal synchrotron 
emission produced by relativistic blast waves which decelerate and radiate by sweeping up particles from 
the circumburst medium (CBM). The terms in this function are derived from basic blast- wave physics. 
Differences in the emission properties of blast waves due to changes in parameter values are examined in 
Section 3. In particular, the strong-dependence of the blast-wave spectra on r is noted. In Section 4, 
we consider the limitations of and requirements for telescopes to be capable of detecting clean and dirty 
fireballs in view of the emission properties of blast waves with different baryon loading. The predictions 
of clean and dirty fireballs as new classes of astrophysical transients are discussed in Section 5, where the 
central results of this study are summarized through a simple rule which governs the overall emission 
properties of blast waves produced by relativistic fireballs. 

2. Spectral Power and Spectral Power Flux of GRBs in the Blast Wave Model 

GRB emissions are most conveniently described in terms of their measured vF v spectral power fluxes. 
For simplicity, we consider only uncollimated, spherically expanding blast waves here. The spectral power 
P(e,t) = vL v is related to the spectral power flux through the relation P(e,t) = A-Kd 2 L vF u , where dh is 
the luminosity distance. Due to the strong Doppler beaming when r > 1, most of the observed flux is 
produced by emission regions located within an angle 9 < l/T of the line-of-sight direction. As the blast 
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wave decelerates into and through the nonrelativistic regime, the observer receives additional contributions 
from emission regions at larger angles to the line-of-sight (Panaitescu & Meszaros 1998a; Waxman 1997b; 
Sari 1998; Chiang & Dermer 1998). This introduces corrections to the above relation between the produced 
power and observed flux, but this effect is not important when r>l, and is neglected here. 

In this limit, we model the temporally-evolving GRB spectrum with the function 

ergs 



P(e,t) 



e x 



s — e 



(l + v/8)P p (t) (i) 



[e/e p (t)]-v + (v/5)[e/e p (t)} 

where, on a vL v plot, the rising slope has index v (upsilon), and the descending slope has index 5 (Sown). 
Both v and S are > 0. The dimcnsionless photon energy e is in units of m e , where particle masses are in 
energy units. Obviously P(e p , t) = P p (t). Expression (1) applies only to the nonthermal synchrotron portion 
of the spectrum, in accord with detailed treatments (Chiang & Dermer 1998; Panaitescu & Meszaros 
1998b) which show that the SSC contribution to the total power output is generally much weaker than 
the nonthermal electron synchrotron contribution. The SSC component can still represent the dominant 
emission process in certain frequency ranges, for example, at 3> 100 MeV gamma-ray energies or at soft 
gamma-ray energies in the afterglow phase when the high-energy flux is at a very low level. At MeV 
energies and below, however, the SSC component can usually be neglected. Compton processes involving 
external photons can also be shown to be negligible unless the surrounding radiation fields have extremely 
large energy densities ;> 10 4 eV cm~ 3 (see, e.g., Dermer, Sturner, & Schlickeiser 1997). 

We assume that the density of the CBM surrounding the location of the burst event can be 
parameterized by a power law in distance from the center of the explosion. Thus 

n(x) = nox- 1 ' . (2) 

The dimensionless spatial coordinate x represents the distance from the explosion site in units of the 
deceleration radius 

= I?^] 173 . (3) 

which characterizes the distance at which the blast wave has swept up « 1/Tq times the initial baryonic 
mass of the fireball (Rees & Meszaros 1992; Meszaros & Rees 1993b). 

It is useful to parameterize the evolution of the bulk Lorentz factor of the blast wave by the expression 

(T , if < ar < 1; 
T(x) = (4) 

[r x-s , if i <x < rl /9 . 

The index g refers to the deceleration regime. If very little of the swept-up energy is radiated over the time 
scale on which the blast wave decelerates, then the blast wave evolves in the non-radiative (or adiabatic 
regime) with g -> 3/2 when rj = (Blandford & McKee 1976; Paczynski & Rhoads 1993). This has a 
self-similar solution independent of To when x » 1. If , on the other hand, the bulk of the swept-up energy 
is promptly radiated away, then the blast wave evolves in the radiative regime with g — > 3 when rj = 
(e.g., Meszaros, Rees, & Wijers 1997 and references therein). Regimes intermediate to these two limits have 
been treated analytically (Dermer & Chiang 1998) and numerically (Chiang & Dermer 1998; Panaitescu & 
Meszaros 1998c). 

Corresponding to Xd is the deceleration time scale in the observer's frame, given by 



t d = 



(1 + z)x d _ (1 + z) 
I> cT 8 Q /3 



(3 - y)E 
4TTnoirip 



1/3 



(5) 
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Eqs. (3) and (5) hold when r\ < 3, and are readily generalized for more complicated radial structure 
involving, for example, inner and outer boundaries of the CBM, or more general radial distributions with 
77 = r](r). In the most general case, rj — rj(r). A more complicated radial dependence of rj than given through 
cq. (2) would produce additional time-dependent variations in the spectral output of the decelerating blast 
wave than derived here. For example, short time scale variability in GRB light curves may result from 
clouds surrounding the GRB source (Dermer & Mitman 1998). The dimensionless spatial coordinate x is 
related to the observed time t through the expression 

' t/t d , if < t < t d and < x < 1; 

X=\ /-IN ( 6 ) 

[(2 ff +l)^-2 3 ] 1 /(2 9 +D , Ht d <t< ^ [Tl +a +2g) and 1 < x < rj /s . 



The existence of a low-energy cutoff in the electron distribution function is reasonable, given that the 
blast wave sweeps up electrons and protons from the CBM with Lorentz factor T in the comoving frame. If 
an efficient mechanism transfers energy from the protons to the electrons, then the low-energy cutoff of the 
electron Lorentz factors can reach values as large as 

7c,min = ie{m p /m e )T (7) 

(e.g., Meszaros et al. 1994). The term £ e < 1 is an electron equipartition factor. Synchrotron radiation from 
an electron distribution with a low-energy cutoff produces a low-energy spectrum with v = 4/3 (Katz 1994; 
Tavani 1996; Cohen ct al. 1997). 

Fermi processes in the blast wave shock can additionally accelerate a nonthermal component of 
electrons and protons. If the injection index is steeper than 3, or if the injection index is steeper than 2 and 
cooling processes arc rapid, then the vF v spectrum falls with 6 > 0. It is assumed here that the values of v 
and (5 and, by implication, the spectrum of the underlying particle distributions do not evolve with time. 
This is clearly not the case due to particle cooling through radiative processes (see Chiang & Dermer 1998 
for a numerical treament). These processes must be carefully considered for interpretations of and fitting 
to afterglow spectra. But these details are not so important for the overall energetics of the blast-wave 
emissions considered here. 

Another crucial uncertainty of blast-wave models is the magnetic field strength H in the comoving blast 
wave frame. This is commonly given in terms of an "equipartition field" obtained by equating the magnetic 
field energy density with the nonthermal particle energy density of the swept-up particles downstream of 
the forward shock. Thus H is parameterized by the expression 

H(G) = ^Trmpnix^Hir/^rix) = 116 [n(x)^(r/4)] 1 / 2 r 300 , (8) 

where n(x) is the density of the swept-up gas in the stationary frame of the GRB source given by eq. 
(2), r is the compression ratio, T3oo = r(x)/300, and is the magnetic-field equipartition parameter. 
Field generation through dynamo processes could strengthen H with time (see, e.g., Ryu & Vishniac 1991; 
Meszaros, Laguna, & Rees 1993 and references therein); flux- freezing and reconnection could weaken it 
(e.g., Meszaros et al. 1994). The choice £ H ~ 1 does not produce GRB spectral forms similar to those 
observed, as excessive cooling produces a spectral component with v = 0.5 (Sari, Piran, & Narayan 1998; 
Chiang & Dermer 1998). Values of £ H <J 10~ 4 are required to give good fits to GRB spectra during the 
prompt phase. Although there is no justification for treating as time-independent when ^ 1, we do 
so here for simplicity. 
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In a tangled field with mean field strength H given by eq. (8), the low-energy electron Lorentz- factor 
cutoff (7) gives a peak in the spectral power at energy H^j 2 min / H cr in the comoving frame, where 
H cr — 4.414 x 10 13 G. To the observer, this energy is boosted and rcdshiftcd by a factor w T/(l + z). Thus 
the observed peak energy e p at time t is given by 

( x-il 2 , if < x < 1; 
e p {t)=£ Q [Y(x)/T fx-^ 2 =£ Q \ (9) 

I x -4 9 -r,/2 ) if i < a; < T y\ 

where 

3.0 x lQ-*nl /2 qTt 
£ ° ~ (lT7) ' (10) 

and 

9 = [Z H {r/A)]^e c . (11) 

The coefficient in eq. (10) is (m p /m e ) 2 (32Trm p ) 1 ^ 2 /H cr = 2.97 x 10 -8 , and So represents the observed 
photon energy of the peak of the vF v spectrum when the blast wave passes through the deceleration radius 
at x = 1. 

The power in swept-up particle kinetic energy in the comoving frame is 

E kc = m p cB(x) F(x) [T(x) - 1] n(x) A(x) (12) 

(Blandford & McKee 1976), where A(x) = A a x 2 = (Airx^x 2 is the area of the spherically expanding blast 
wave and B(x)c = [1 - r(x)- 2 ] 1/2 c is its speed. If £ is defined as the fraction of swept-up energy retained 
in the comoving frame, then (1 — Q is the fraction that is dissipated by the blast wave shock. Although this 
dissipation could be in the form of energetic particles which escape directly from the shock, or in the form 
of bulk kinetic energy which is randomized and stored in the reverse shock to be radiated away later or 
used to accelerate the blast wave front (see, e.g., Meszaros & Rees 1997; Panaitescu & Meszaros 1998a, c), 
we assume that essentially all of the dissipation occurs in the form of photon radiation. 

As previously discussed, only the emitting region along the line of sight makes a significant contribution 
to the observed flux in the limit r > 1. The observed power is boosted and redshiftcd by the factor 
T 2 /(l + z) 2 over the radiant power in the comoving frame. Consequently we obtain an expression relating 
the observed power to the swept-up power, given by 

^rad - n r ' , 2 (1 - C)Ac - (1 - C)m p cT' l n(x)A(x)/(l + z) 2 . (13) 

(l + zy 

The spectral power given by eq. (1) satisfies the normalization 

poo poo poo 

E = dtE lad (t)= dt dee- 1 P(e,i), (14) 
Jo Jo Jo 

where E is the total energy in baryons at the end of fireball coasting phase. Thus the radiated power 

E Tad = +6- 1 )P p (t) . (15) 

We use the approach of Dermer & Chiang (1998) to find that in the regime 1 <i< Cr^ 9 , the fraction 
of swept-up energy that is dissipated promptly is given by (1 — () = (2g — 3 + r])/g (see Appendix), and we 
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assume that this relation holds elsewhere as well. Consequently we derive an expression for the measured 
bolometric spectral power, given by 

(x 2 ^, if < x < 1 ; 

P P (t) = n J (16) 

[ x 2- n -4g^ < x< rl /9 . 

The coefficient 

_ (2g-3 + ^myd^no A) c(2g - 3 + ??)rg /3 1/3 , 2/3 2 / 3 , . 

Ho - 2 fl («-i+«-i)(l + z)» " 2 fl («-i + J-i)(l + z)* (47rmp " 0) (3 _ V) E ° 

represents the bolometric luminosity (ergs s _1 ) at the time td of peak power output. We again note that a 
more complicated expression would hold if r\ were not constant, thereby rendering £ dependent upon time. 

It is straightforward to determine the time-dependence of e p (t) and P p (t) in different regimes from eqs. 
(6), (9), and (16). For the parameterized density distribution (2), the peak photon energy e p (t) oc t~ v ^ 2 
when t < td, whereas for t > td, e p (t) oc t _ ( 12 + ? ')/ 8 and oc /j-( 24 +*?)/ 14 in the non-radiative and radiative 
regimes, respectively. For a uniform medium (77 = 0), e p (t) is constant for t < td, whereas for t > td, 
e p (t) oc t~ 3 / 2 and oc t~ 12 / 7 in the non-radiative and radiative regimes, respectively. The peak spectral 
power P p (t) increases oc t 2 ~ v for t < td, and P p (t) oc t~ 1 ~ r >/ 4 and oc t _ ( 10 + 1 ')/ 7 in the limiting non-radiative 
(g = 3/2) and radiative (g = 3) regimes, respectively, when t ^> td- For a uniform medium with rj = 0, 
P p (t) oc t 2 for t < td, and P p (t) oc t -1 in the non-radiative regime when t > td (Sari 1997). In the 
radiative regime, P p (t) oc i~ 10 / 7 when t > td- From these relations, one can see that the power emitted 
in lower-energy radiation such as the optical band will always be less than in the higher energy X-ray and 
gamma-ray bands for declining density profiles (i.e., i] > 0). 

One can also easily determine the temporal indices measured at a fixed observing energy e by consulting 
equations (1), (6), (9), and (16). Let t p (e) represent the time when the spectral power flux at e reaches its 
maximum value. It is found by solving e = e p (t p ) using eqs. (6) and (9), yielding 

2a+i 

(18) 



«p(e)= !i 



2.9+1 



" 4g + r,/2 

+ 2g 



Note that t p (e) — ► td when e 3> £0, whereas t p (e) ^> t p when e <C £q. The displacement of the pulse peaks at 
lower energies may act alone or in concert with synchrotron cooling (Chiang 1998; Dermer 1998) to produce 
the observed energy dependence of the GRB pulse width (Fenimore et al. 1995; Piro et al. 1998b). Effect 
(18) represents a prediction of this model, and will be explored in more detail elsewhere (Dermer, Bottcher, 
& Chiang 1998). 

Two regimes characterize the time-dependence of the flux at high energies e > £q, namely 



e 



{t/td) 2 -^ 8 ' 2 , if *«*«!, 



PM) = (i + -)n (-)- d x <^ (19) 

1 [(l + 2. 9 )t/t d ][ 2 -"(l + *72)-4 3 (l+5)]/(2 3 +l) j iH>7 ^ 

In contrast, three regimes characterize the time-dependence of the flux at energies e <C So- They are 

f (1 + v/s)(f o r(t/t d ) 2 -^ 2 , if t « t d , 



P(e,t) = U x < 



(1 + v/S)(^) v [(l + 2 5 )t/t d ][ 2 -'Ki-W2)-4 3 (i-.)]/(2 3 +i )j if u ^ t <<; y e); (2Q) 
_ (1 + + 2g)t/tdf-^ 1+s ^- 49(1+s ^^ +1 \ if t > t p (e). 
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Note that the late time behaviors in eqs. (19) and (20) are identical. 

At times t > td, two types of behavior occur depending on whether the observations take place at 
e <C £o or e ^> £ . In the former case, there are two branches in the time profile related to the time t p (e) 
when e p sweeps through the observing energy e. When td -C t <C t p (e), 

P(e,t)ozt X1 where X i = [2 - - v/2) - 4g(l - v)]/(2g + 1) . (21) 

When t S> t p (e), the time profile decays according to the relation 

P(e,t)cct-^ where X i = [4.9(1 + $) + r,{\ + 5/2) - 2]/{2g + 1) (22) 

In the latter high-energy case (e 3> £o), a single power-law decay in the time profile is found for t ^> td 
which follows the behavior given by eq. (22) . Table 1 lists the derived temporal slopes \i an d X2 for various 
values of v, 8, and g. 



3. Model Spectra 

Table 2 gives a list of the parameters which, when placed in eqs. (5), (6), (9), (10), (16), and (17), 
can be used to calculate temporally-evolving model spectra using spectral form (1). Fig. la shows such a 
calculation employing the parameter values listed in Table 2 and called "standard" henceforth. For these 
parameters, the deceleration time scale td = 58 s and £q = 0.4, giving a duration and a 200 keV vF v peak 
energy early in the event which are characteristic of typical GRBs. The smoothly joined broken power-law 
form resembles the canonical GRB "Band" -type spectrum (Band et al. 1993). This model burst radiates a 
peak power n = 3.8 x 10 50 ergs s^ 1 at t d and would be below the BATSE flux threshold at redshift z = 5, 
but would be strongly detected if the GRB source were at z = 1. 

The spectra in Fig. la are plotted at observing times from 1 /is to 10 7 s in factor-of-10 increments. 
The solid lines which overplot the spectra connect the peak energies e p (t) and peak luminosities P p (t) 
at different observing times. When t < td, e p (t) is constant because 77 = in this example, whereas it 
declines at later times according to the relation e p (t) oc /j _4 s/( 2 3+ 1 ) oc i~ 1,6 , noting eqs. (6) and (9) with 
g = 2. For these same parameters, the peak power output P p (t) oc t 2 when t < td, but P p (t) decreases 
cx f(2-4g)/(2 9 +i) K £-6/5 wnen i ^> ^ noting eqs. (6) and (16). We see that observations of the temporal 
decay of the vF v peak relate the indices g and 77 which characterize the radiative regime and density profile 
of the external medium, respectively. Again, the possible spatial dependence of r\ complicates any simple 
interpretation of the measured behavior. 

Fig. lb is a representation of the model data shown in Fig. la in the form of time profiles at different 
observing energies. Light curves are shown ate = 10~ 12 , 10~ 9 , 10~ 6 , 10~ 3 , 1, 10 3 , and 10 6 , roughly 
corresponding to radio, mm, optical, X-ray, MeV, GeV, and TeV frequencies, respectively. The light curves 
at t > td display one or two power-law behaviors depending on whether e > £0 = 0.4 or e < 0.4, respectively, 
for the reasons given in the previous section. For the very late time behavior with t 3> t p [e) > td, the 
temporal decay index \2 = [4<?(1 + S) — 2]/(2g + 1) = 1.52, in agreement with the behavior shown in Fig. 
lb. It is interesting to note that temporal indices in the observed range 1.1 <; \ ^ 1-5 require flat spectra 
with 8 £ 1/3 when rj = (see Table 1) . 

The GRB spectral indices v and 5 are assumed for simplicity not to evolve with time. Thus spectra at 
photon energies e -C e p and e » e p rise and fall with the indices v = 4/3 and 8 = 0.2, respectively, for the 
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parameters given in Table 2. To examine the effects of different parameter values on the observed spectra, 
it is therefore sufficient to plot the P p (t)-e p (t) trajectory at various times, as indicated in Fig. la by the 
straight lines. Figs. 2a-e presents a parameter study in this form. The Table 2 standard is shown in each of 
these figures by the solid lines with open circles at times incremented by factors-of-10 starting from 1 us. 

Fig. 2a shows the effects of changing the index which characterizes the radiative regime. Here we 
show two cases with g = 1.6 and g = 2.9, in addition to the g = 2 standard. The case g = 1.6 is near the 
non-radiative limit (if g = 3/2, there is no radiation by definition), and g = 2.9 is near the fully radiative 
limit. The weakly radiative blast wave with g — 1.6 persists for a longer period of time, though it radiates 
less power in its early phases. The highly radiative blast wave with g — 2.9 is nearly an order-of-magnitude 
more luminous in its early phases than the g = 1.6 case, but its inertia drops more quickly and it therefore 
has a shorter luminous life than in the less radiative cases. 

At this point, a remark concerning a serious issue regarding lack of self-consistency of the parameter list 
is in order. The type of radiative regime specified by g essentially gives the fraction of swept-up power that 
is dissipated. If nonthermal synchrotron radiation represents the dominant power drain from the blast wave, 
as we assume, then a fraction « 2 — g~ (3 — 77) of the incoming swept-up power is promptly transformed into 
radiation. This places demands on the underlying particle distribution and comoving energy densities which 
are not considered here, but must be kept in mind in any analytic treatment employing a parameterization 
of blast-wave dynamics according to an expression like eq. (4). A detailed numerical simulation with a 
self-consistent treatment of plasmoid/blast-wave dynamics is essential to model spectra accurately (Chiang 
& Dermer 1998). 

Figs. 2b and 2c show the effects of varying the density n and the explosion energy Eq, respectively, over 
four orders of magnitude. In higher density environments, a blast wave's energy is radiated more rapidly 
and, provided q is the same, in the form of higher energy photons; thus its peak power output obtains larger 
values than in more dilute environments. Fireballs with larger total energies but with equal baryon-loading 
persist for longer periods of time because both td and the time scale to decelerate to nonrclativistic energies 

1 /3 

are oc E ' (see eqs. [5] and [6]). For a fixed baryon loading factor (or constant r ), the overall duration of 
the GRB is multiplied by a factor oc (Eq/uq) 1 ^ 3 as a consequence of the deceleration time scale (see eqs. [3] 
and [5]). This dependence comes from dimensional analysis in scaling solutions for spherical hydrodynamics, 
and is present in both the Sedov solution (see, e.g., Lozinskaya 1992) and its relativistic equivalent derived 
by Blandford & McKee 1976). The peak power output changes to compensate for the different duration; 
thus P p (td) oc E^J 3 when Eq varies and no is held constant, and P p {td) oc n^ 3 when uq varies and Eq is 
constant (eq. [17]). This latter dependence is a function of the radiative efficiency, which is assumed to be 
identical to the factor (1 — £) derived in the Appendix. Consequently, the dimensional dependence of P p (t) 
is more complicated than that of td- In any case, the durations and powers are not strongly dependent on 
changes of n and Eq . 

Fig. 2d shows the effect of varying the equipartition term q over four orders-of-magnitudc, from 
q = 10 -5 to q = 10 _1 . A change in q by a constant factor causes the curves to be displaced horizontally by 
the same constant factor, since c p (t) oc q but P p (t) is independent of q. 

The effects of different amounts of baryon loading are shown in Fig. 2e. The spectra are found to be 
extremely sensitive to changes in r . The mean photon energy Eq = e p {td) at the time of peak-power output 
Ho is oc Tq. Thus most of the radiation of a clean fireball with Tq = 3000 is in the form of photons with 
energies eight orders of magnitude larger than the emission from a dirty fireball with To = 30, provided 
that q remains constant. The duration of peak power output is just the deceleration timescale td which is 
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g /3 8/3 

oc r o .To compensate for the constant total energy in the fireball, IIo oc T . Thus the duration of peak 
power output of the clean fireball case is a factor ~ 2 x 10 5 shorter than that of the dirty fireball case in 
this example, but the magnitude of the power emitted during this period is a factor of ~ 2 x 10 5 greater. 

Clean fireballs produce extremely short pulses of radiation which are carried by very high-energy 
photons. For example, when T = 3000, the peak power output of « 10 53 ergs is carried by ^ 1 GeV 
photons in a burst lasting <; 0.1 s. Again, the total energy remains roughly constant, though it is carried by 
many fewer, though much higher-energy photons. Its late-time behavior approaches that of fireballs with 
smaller values of To, but with lower power to compensate for the energy radiated during the early phases of 
the event. For the dirty fireball with r = 30, the peak power output at UV energies lasts for ~ 10 minutes 
and develops into a broadband radiation pulse which evolves to lower energies at later times. Incidentially, 
this behavior is found in blazar flares (e.g., Marscher & Gear 1985), though a closer examination must be 
made to see if this model provides an adequate explanation for such phenomena. 

The GRBs which trigger BATSE display a very broad bimodal duration distribution ranging from ~ 30 
ms to « 500 seconds (Kouveliotou et al. 1993; Belli 1995). A short, hard component of GRBs is seen with 
durations ranging from w 30 ms to 1 second, with a mean duration of w 0.2 seconds. A longer duration 
component ranges from w 2-500 seconds, with a mean duration of 15 - 30 seconds. This hardness-duration 
correlation is in accord with the behavior outlined above for blast waves with different baryon-loading, 
where the class of short, hard bursts are produced by fireballs which are less baryon-loaded than average. 
A model for GRBs in terms of a distribution of To is needed, however, to establish whether these GRBs 
constitute a separate class, or an extension of the typical loaded fireball in the cleaner limit (Dermer et 
al. 1998). Different levels of baryon-loading might be expected if GRBs originate from multiple sources, 
as indicated if the claimed association of GRB 980425 with SN 1998bw is correct (Kulkarni et al. 1998b). 
For the purposes of this paper, we consider a continuous range of fireball-loading parameters rather than 
different levels of baryon loading. 

Thus we find that the observed spectra from fireball/blast waves, which decelerate and are energized 
by the process of sweeping up material from the CBM, are most strongly dependent on the baryon-loading 
parameter r and the equipartition parameter q. The term q determines the mean photon energy of the 
radiant emission through the relation £$ oc qT^ (cqs. [9]- [11]). In the most general case, q = q(To), and the 
observational properties of the fireball would depend on the exact form of this dependence. We assume 
that q is only weakly dependent on r , that is, we assume that q depends on r though a function which 
varies much less strongly than oc T^ 4 . In what follows, we examine the observational signatures of fireballs 
with different baryon loading r but with constant q, keeping in mind the possibility of an underlying 
dependence of g on Tq. 



4. Prospects for Detecting Clean and Dirty Fireballs 

The total number of photons collected from a GRB with energies near the mean photon energy of peak 
power output is roughly given by 

&TTd z L m e £o 

where A<j et is the effective area of the detector at energy £q. The deceleration time scale (5) defines the 
rough duration of peak power output, and eq. (10) gives the mean photon energy of the peak power output 
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at t = td- Hence 

2/3 



t d (s) = 1730 (j-^j [(1 - r?/3)(l + z)E 5i ] 1 / 3 - 31 (|^) 2 / 3 [(1 - r 7 /3)C 6 ^54] 1/3 



(24) 



where g_ 3 = g/10 3 and ( & = (z + l)/6. Recalling the expression for the magnitude of peak-power output 
given by eq. (17), we obtain the result 

_ (2. 9 -3 + 7?) (3-r?)£ A det _ E 54 A det 
7 ~ 2g(v-i + + «) 47rm e d 2 £ 7 4^(1 + ^) ' 1 J 

where c?29 = ciL/(10 29 cm). Note that c/29 = 1 when z = 5 in a qo = 1/2 cosmology with a Hubble constant 
equal to 65 km s _1 Mpc -1 . For our standard parameters r\ = 0, v = 4/3, and S — 0.2, 

r 0.32, for .9 = 1.6; 
K 7 = I (26) 
[ 2.4, for .9 = 2.9. 



4a. Clean (r > 300) Fireballs 

The 100 MeV threshold of EGRET corresponds to £ ~ 200. To produce 100 McV nonthermal 
synchrotron photons mostly from electrons with energies ~ m e "f e _ m i n requires fireballs with baryon-loading 

T = 76[^-^l}^ = 1420(^— )V< , (27) 
n ' q n 2 ' 9_ 3 

where we set So = 200 on the right-hand-side of eq. (27). The number of photons detected from a clean 
fireball, given the standard parameters used in eqs. (25) and (26), is iV 7 rs £? 5 4Ad e t[cm 2 ]/[3750d29C6] and 
iV 7 « _E 5 4 J 4det[cm 2 ]/[490c?29C6] for 9 = 1-6 and g = 2.9, respectively. To detect clean fireballs requires that 
7V 7 > 1. Thus sources of clean fireballs must be located at distances d 2 Q given through the relation 

r 0.5 for .9 = 1.6 

Cl' 2 d 2 9 <; i^iooo <^ , (28) 

[ 1.4, for .9 = 2.9 

where 1000+iooo cm 2 is the effective collecting area of the detector at £q s=s 200. 

The effective area of EGRET at 100 McV in the wide field mode (e.g., Kurfcss et al. 1997) is - 800 
cm 2 on-axis and ~ 400 cm 2 at 20° off-axis, so clean fireballs would have to be located at z -C 5 and at 
z ^ 5 to be detected with EGRET in the non-radiative and radiative limits, respectively. Bursts from 
clean fireballs with £ ~ 200 would have durations t d w 0.9[9 2 _ 3 Ce(l — rj/fyE^] 1 / 3 (eq. [24]) at z = 5 and 
shorter durations when z -C 5 and 9_3 <C 1. The spark chamber of EGRET requires « 0.1 s to recover 
from a photon event, although it can detect multiple photons if they arrive within a few microseconds. 
Consequently, EGRET is insensitive to more than one photon arriving between several microseconds and 
~ 100 ms, which is a critical range for the clean fireballs considered here. 

It might be surprising that no sub-second bursts of gamma rays have been reported (Fichtel et al. 
1994; Fichtel & Sreekumar 1997) from some members of the clean fireball class which are located at z <J 5. 
We propose three possibilities to explain why this class of clean fireballs with peak power output near 100 
MeV has not been discovered: 
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• There are many fewer clean fireballs, as defined by eq. (27), than those which produce the loaded 
BATSE GRBs. 

• The duration of clean fireballs is typically <; 0.1 s, so that deadtime limitations have prevented their 



• Members of the clean fireball class evolve through cosmic time differently than the loaded fireballs 
which produce the GRBs detected with BATSE, and they do so in such a way to have escaped 



If the third point does not apply, then we can derive a requirement on the rcdshift of the dimmest 
BATSE bursts to agree with the failure to detect clean fireballs with EGRET. Very crudely, we scale the 
number of clean fireballs which could be detected with EGRET to the BATSE detection rate (~ 800 yr _1 
full sky) and the EGRET ficld-of-vicw (~ l/25th of the full sky) and lifetime (« 5 yrs). If clean fireballs 
can only be detected by EGRET within z <; ze as given by eq. (28), then for EGRET not to have detected 
at least one member of this class, we find that 



The term k c f is a class enhancement factor correcting for the number of clean fireballs relative to the 
number of loaded fireballs which produce GRBs detectable with BATSE. The exponent j represents a 
cosmological scaling, and j = 3 for Euclidean space. Evolutionary and cosmic expansion effects will 

— 1/3 

introduce modifications to j not considered here. Simply taking j — 3, we find that ze <^ 0.2zsfc c y 
for the clean fireballs not to have been detected with EGRET. Noting eq. (28), this implies that even if 
k c f ~ 1, sensitivity and FoV limitations could have prevented their discovery. Adding to the difficulty of 
the discovery is the deadtime limitations given by the second point above. 

Stronger conclusions require a detailed size-distribution study involving the cosmological evolution 
of loaded GRB fireballs and clean fireballs in different scenarios, for example, the hypernova scenario 
(Paczynski 1998; Wijers et al. 1998). Given the uncertainty in the parameters, it remains an open question 
whether the larger effective area of the proposed GLAST instrument will be adequate for the discovery 
of the clean fireball class, or if all that is required is a wide FoV, large area gamma-ray calorimeter with 
fast timing. The discovery of subsecond bursts of ~ 100 MeV photons from clean fireballs remains a real 
possibility which, should such an event be detected, is not to be confused with the Hawking radiation from 
evaporating mini-black holes by virtue of its afterglow behavior. 

Clean fireballs would be difficult to detect with the Solar Maximum Mission, which accumulated on 
time scales of 16 s and was not sensitive to short bursts with durations ^ 2 s (Harris & Share 1998). 
The short hard clean fireball bursts would also be difficult to detect with COMPTEL on CGRO, noting 
its smaller effective area (~ 10 2 cm 2 ; see, e.g., Hanlon et al. 1995). Nonetheless, further study of the 
COMPTEL sample of GRBs detected with BATSE, and a search for unidentified MeV transients will 
provide additional constraints on the number and redshift distribution of clean (r 3> 300) fireballs. 



discovery. 



detection. 




(29) 
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4b. Dirty (r < 300) Fireballs 

Because a large number of X-ray telescopes operate in the 1-10 keV range, we can use £ ~ 3/511 to 
define the dirty fireball regime, which includes those blast waves with 

T ~ 21 I^^] 1 / 4 = 104 (^|M 1/4 (30) 

(compare eq. [27]). Detection of dirty fireballs with X-ray telescopes is difficult because of the intense 
backgrounds. A signal of iV 7 counts given by eq. (23) is detected above a background of Bq counts at the 
n a significance level provided that 

N 7 % n^Ba) 1 ' 2 . (31) 

The number of background counts 

B G S Ae dct ■ $ dif (n, e d ct) ■ Ml ■ At ■ A dct , (32) 

where, from right to left, the terms are the detector area, the observing time, the field-of-view (FoV) of the 
detector, the diffuse photon flux weighted by an average photon energy within the band, and the bandwidth. 
The background noise is composed of instrumental, magnetospheric, zodiacal, Galactic, and extragalactic 
components, which are time-dependent in general. Here we give a background estimate for X-ray telescopes 
considering only the diffuse cosmic X-ray noise background, yielding a lower limit to Be- 
lli the 3-100 keV range, Boldt (1987) gives the expression 

e$XRB(£ ' n) [ ^V7i ] ~ 5 ' 6 [ f|ely r0 ' 29 cx P[-£( kcV V 40 kcV l ( 33 ) 

for the diffuse extragalactic X-ray background flux. At 3 keV, corresponding roughly to the mean energy of 
photons detected by an X-ray telescope sensitive in the 1-10 keV band, the diffuse flux is ~ 5 ph cm~ 2 s _1 
sr _1 . 

Substituting eqs. (23) and (32) into eq. (31) with the stated diffuse flux gives the result that detection 
at the n a level requires an X-ray telescope with 

cm 2 ! 3.4 x 10~ 4 4, n 2 (1 + z) 2 At 

~\ * K 2 El f(At/t d ) ( } 



Aict 

An 



The function f(y) weights the significance of detection by the observing interval over the period in which 
the fireball's emissions are monitored. Its form depends on the time-dependence of the light curve of the 
fireball emissions, which can be modeled using the eq. (1), and the epoch of observation. In an illustrative 
regime, f(y) = y 2 for y < 1, and f(y) = 1 for y > 1. The detection efficiency is clearly optimized when 
At ~ td- The deceleration time scale for a burst peaking at £q — 3/511 = 0.006 is 

t d (s) ~ 970 M_ 3 (l - r//3)£ 5 4] 1/3 • (35) 

The seach for dirty fireballs should be conducted in the range of time scales from 10 2 -10 4 s, recognizing the 
difficulty due to orbital variations of the background radiation environment. 

When r] = 0, one obtains the requirement for an X-ray telescope to be able to detect a dirty fireball at 
the n a level, namely 

A dct . Cl /3 4 9 (nj5) 2 q % 3 j 94 °' «* g = 1.6; 

~mT * W 3 (36) 

^54 [ 125, for g = 2.9. 
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Typical X-ray detectors have effective areas ~ 10 2 cm 2 (see Table 3 for general specifications of some X-ray 
telescopes used to detect GRBs), so we scale to a detector of area 100j4ioo cm 2 . To detect dirty fireballs 
above the diffuse X-ray background means that the detector's FoV must satisfy the relation 



Thus an X-ray detector has to have a FoV <; 10° x 10° or <; 80° x 80° in order to detect dirty fireballs at 
z ~ 5, depending on whether the blast wave decelerates primarily in the non-radiative or radiative regimes, 
respectively. Fireballs at z C 1 place much weaker constraints on the detector's FoV, as indicated in 
equation (37), and then it becomes a question of the number of such sources one could expect to see. 

Small FoV instruments will only detect an X-ray flash from a dirty fireball very rarely unless the class 
of dirty fireballs is populous. Given that there are ~ 800 BATSE-detectable GRBs per year averaged over 
the full sky, this implies a chance probability of ~ 0.02 • kjf ■ k\ per year for a detector with a k° d x k d FoV. 
The class enhancement factor kdf of dirty fireballs compared to loaded BATSE/GRB fireballs is not known 
at present, though it can be strongly constrained by modeling the duration distribution of GRBs detected 
with BATSE or with other burst detectors (Dermer et al. 1998). Searches for dirty fireballs should be made 
with the Wide Field Camera data on Beppo-SAX, and with future wide FoV imaging X-ray detectors. 
Concrete deductions about the probability of detecting serendipitous dirty fireball X-ray flashes requires a 
size distribution study tied to a particular detector's characteristics. 

It is feasible that some of the X-ray flashes observed in archival searches of HEAO (Connors, 
Serlemitsos, & Swank 1986) and ROSAT PSPC (Li et al. 1998; Sun et al. 1998) data are emissions from 
dirty fireballs only if k df > 1, given that the effective FoV of the ROSAT PSPC is - 1° x 1°. The GRBs 
detected with Ginga (Strohmayer et al. 1998) have vF v peaks which are on average at lower energies than 
those detected with BATSE. Because Ginga was sensitive to lower energies events than BATSE, this might 
indicate that it is more sensitive to a population of GRBs with greater baryon-loading. A straightforward 
prediction from this model is that the mean duration of detected GRBs is inversely correlated with the 
triggering energy range. 



A slewing strategy such as the one so successfully demonstrated by the Beppo-SAX team and 
collaborators (e.g. Costa et al. 1997; van Paradijs et al. 1997; Piro et al. 1998a,b) opens rich possibilities for 
discoveries which must be carefully weighed in GRB telescope design. Here we have attempted to provide 
an exposition of the blast-wave model which provides a simple way to model the evolving spectral power 
fluxes, and which furthermore helps us to understand dedicated and serendipitous GRB observations. 

Eq. (1) is the proposed time-dependent spectral form. It employs the 9 parameters listed in Table 
2 that follow from the basic blast-wave model. Its crucial underlying assumptions are that the photon 
indices v and 5 and the equipartition term q are constant in time, and that q is independent or only weakly 
dependent on r . The energies of photons which carry the bulk of the power are oc gFg and are therefore 
most strongly dependent on the value of Fo- Over this looms a question of lack of self-consistency, since 
specification of the radiative regime g constrains the magnetic field and photon energy densities and particle 
distributions to produce the requisite radiative power. Comparison of eq. (1) with the results of a numerical 
code employing a self-consistent treatment of plasmoid dynamics (Chiang & Dermer 1998) indicates that 




for g = 2.9. 



for g = 1.6; 



(37) 



5. Discussion 
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the proposed analytic description of the blast wave emission spectra remains accurate except when the blast 
wave is highly radiative. 

In spite of these concerns, eq. (1) can be easily employed to make estimations of imaging capabilities 
of burst detectors given slewing rates of narrow field instruments toward a cosmic transient. It can be used 
to calculate photoelectric absorption variations in the GRB environment and its imprint upon the observed 
burst spectrum (Bottcher et al. 1998). It can also, of course, be used to model the temporal behavior of 
GRB spectra and to examine GRB phenomenology in terms of the fireball/blast- wave model (Dermer et al. 
1998). 

By displaying the spectral power flux in the form of eq. (1), we obtained an important rule on fireball 
emissions relating the observed peak power n , the characteristic duration td of this peak power output, and 
the mean energy £q of photons which carry this energy. The characteristic time scale over which the bulk 
of the power is radiated is the deceleration time scale td (eq. [5]; Rees & Meszaros 1992), and td oc r o 8//3 , 
where Tq specifies the baryon loading. The peak power output IIo (eq. [17]) is also oc T^ 3 , assuming that 
the bulk of the dissipated energy is in the form of photon radiation. Because most of the total particle 
kinetic energy, E , originally in the fireball during its coasting phase is radiated during the observed time 
td, we have 

E w n t d . (38) 

Precisely, 

W' = *»i^( lisJ , ;(1+ili for9=29 « 39 » 

(In the evaluation on the right-hand-side of this expression, we let r\ = 0, v = 4/3, and S = 0.2.) Moreover, 
the bulk of this energy is observed in the form of photons with characteristic energy £q given by eq. (10); 
note that Eq oc qT^. For these parameters, a full 25% of the energy of the GRB source is radiated within 
the deceleration time scale for a highly radiative fireball. 

Eqs. (38) and (39) can be used to summarize the emission properties of fireballs with different baryon 
loading and prospects for detecting such fireballs. Fig. 3 shows the variation of n , td, and £ when r 
ranges from 100 to 1000. The dependence of these observable due to variations in E 5 4 and n are also 
indicated. The overall powers, time scales, and photon energies of peak power output of fireballs with 
To ~ 300, and for clean and dirty fireablls can be determined from this figure. 

We used benchmark characteristics of X-ray, soft gamma-ray, and medium-energy gamma-ray telescopes 
to establish criteria for clean (r 3> 300), loaded (r ~ 300), and dirty (r <C 300) fireballs. Detection of 
clean fireballs at ;> 100 MeV energies implies r ;> 1000 (cq. [27]), and the detection of dirty fireballs at 
~ 3 keV energies implies r < 100 (eq. [30]). The clean fireball class was not discovered with EGRET; 
this means either that there are relatively few sources of energetic clean fireballs, or that the powers and 
durations of clean fireballs are such that sensitivity and deadtime limitations of previous > 100 MeV 7-ray 
detectors have prevented their detection. 

Here we qualify our use of eq. (38) which gives the general observational properties of fireballs with 
different baryon loading parameters, EGRET has, of course, detected > 100 MeV gamma-ray emission 
from GRBs detectable with BATSE from the loaded (100 £ T £ 1000) fireball class. Many of these sorts 
of GRBs will undoubtedly be detected with the GLAST mission as well. But the > 100 MeV and GeV 
photons detected to date (e.g., Hurley et al. 1994; Dingus 1995) are not those produced by the electrons 
near the low-energy cutoff of the electron distribution function. These photons are much more likely to be 
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the SSC emission (e.g., Chiang & Dermer 1998) or the nonthermal synchrotron radiation from ultra- high 
energy protons (Vietri 1997b; Bottcher & Dermer 1998) accelerated in the loaded fireball blast wave. These 
emissions are distinguished by their slow decay, which explains why EGRET often detects this emission long 
after the main portion of the BATSE burst has decayed. What EGRET has not discovered but is implied 
by the fireball/blast-wave model is a class of clean fireballs which produces a luminous subsecond burst of 
3> 100 MeV radiation. The prediction of such a discovery with improved 7-ray telescopes is made here. 

We also used analyses outlined by eq. (38) to describe the properties of dirty fireballs. Because IIo and 
£0 are lower in dirtier fireballs, the diffuse sky backgrounds make detection of such a class of objects by 
wide-field instruments problematic. On the other hand, the expected rate of detection of dirty fireballs by 
pointed instruments with a few square degree FoV is very unlikely unless the rate of dirty fireball explosions 
greatly exceeds the rate of loaded BATSE/GRB fireballs. Whether any members of this class has been 
discovered is uncertain, though X-ray flashes have indeed been detected in the HE AO and ROSAT data 
bases. We predict that moderate FoV instruments (~ 1% of the full sky) with surveying capabilities will 
discover dirty fireball X-ray transients with typical durations of 10 2 -10 4 seconds. 

Beppo-SAX has shown that improved imaging and rapid response are key to making progress in 
solving the GRB mystery. Granted this success, the design of a new GRB mission should therefore employ 
a hard X-ray/soft GRB detector and a narrow FoV X-ray telescope which can rapidly slew to the GRB 
position. But to discover the sister classes of clean and dirty fireballs implied by the fireball/blast-wave 
model will require additional X-ray and gamma-ray capabilities as outlined above and considered in more 
detail elsewhere (Dermer et al. 1998). 
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Observatory Guest Investigator Program. CD thanks Jim Kurfess, Neil Johnson, Mark Strickman, and 
Gerry Share for enlightening discussions, and the referee for a useful critique. The work of JC was performed 
while he held a National Research Council - Naval Research Laboratory Associateship. MB acknowledges 
support by the German Academic Exchange Service (DA AD). 



Appendix A. The Fraction of Radiated Energy 

An expression is derived in the asymptotic regime x 3> Xd which relates the fraction Q of swept-up 
energy retained in the blast wave to the index g which characterizes the radiative regime and evolution of 
the Lorentz factor T time scale of the blast wave (see eq. [4] ) . We start from the equation of momentum 
conservation (eq. [6] of Dermer & Chiang [1998]) and note that the energy in swept-up particles accumulated 
up to the point x is equal to 



00 X 

J d P7 A p ( 7 ) = J dxn ext (x)A(x) [CT(i) + (1 - C)] ■ 



(Al) 





In the asymptotic regime X4 <C x <C x tr , where x tr = XdT^ 9 and T(x) = r (x/xd)~ 9 3> 1, this yields 
r'(x) g 



n ° A ° fey 



r 2 (£) xT(x) 



N ai + n A 



r oC/^(t) +(W)p(^) 



2-rj 



(A2) 
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Unless the blast wave is in the extreme radiative regime, the denominator on the right-hand side of eq. (A2) 
is dominated by the first integral. Thus if x -C C^tr, eq. (A2) becomes 

*r,2^^y-, ( A3) 



T x d \x d J (Tox d \x d 

or 



C=^^. (A4) 



Table 1 : Slopes of the Temporal Evolution of the Flux at Fixed Photon Energy when t^> t d and r\ = 



slope 


g 


v = 1 


u = 4/3 












Xi 


1.6 


0.48 


0.98 












Xi 


2.9 


0.29 


0.86 
















5 = 


5 = 1/5 


S = 1/3 


S = 1/2 


5 = 3/4 


5 = 1 


5 = 2 


X2 


1.6 


1.05 


1.35 


1.56 


1.81 


2.19 


2.57 


4.09 


X2 


2.9 


1.41 


1.75 


1.98 


2.26 


2.69 


3.12 


4.82 



Table 2: Standard Parameters for Model of Evolving GRB Spectral Power Flux 



Parameter 


Standard Value 


Description 


To 


300 


initial bulk Lorentz factor 


? 


10 -3 


equipartition term 


9 


2 


index of F evolution 


Eq 


10 54 ergs 


total particle kinetic energy in GRB 


V 


4/3 


spectral index of rising portion of vL v spectrum 


8 


0.2 


spectral index of falling portion of vL v spectrum 


z 


5 


cosmological redshift 


no 


10 2 cm" 3 


density at x = Xd 


V 





index of density distribution 



Table 3: Characteristics of Operating Missions with GRB Detectors 





Energy band (keV) 


FoV 


Effective Area (cm 2 ) 


Sensitivity b 


Imag 


ing 


CGRO BATSE 


50- 


-300 


4tt 


1800 cm 2 a 


6 x 10~ 8 c 


few de 


grees 


SAX GRBM 


40- 


- 700 


20° x 20° 


120 




10'- 


>3' 


SAX WFC 


1.5 


- 26 


20° x 20° 


120 


- 10- 10 in 10 3 s d 


10'- 


> 2' 


SAX NFI LECS 


0.1 


- 10 


0.5° 


22@0.25 keV 


- 3 x 10~ 14 in 10 5 s 


50" 


SAX NFI MECS 


0.1 


- 10 


0.5° 


150@6 keV 









"Each Large Area Detector 
6 ergs cm -2 s _1 

"Large Area Detectors; 10 s GRB 

d Beppo-SAX data from homepage; see also Piro et al. (1998) 
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Fig. 1. — Spectra and light curves for a model GRB using the standard parameters in Table 2. (a) Evolution 
of model GRB spectra at different observing times from 1 /is to 10 7 s are shown in factor-of-10 increments. 
Solid lines plot the trajectory of the peak power P p (t) as a function of the photon energy e p (t) of the peak 
of the vL v spectrum at different observing times, (b) GRB light curves at different observing energies for 
the model GRB spectra shown in (a). 




Fig. 2. — Trajectories of model GRBs in the P p (t)-e p (t) plane. All parameters are given in Table 2, except 
as noted, (a) Trajectories for different values of the index g characterizing the radiative regime and evolution 
of T through eq. (4). 




Fig. 2.- 



(b) Same as Fig. 2a, but for different values of the density n of the circumburst medium. 




Fig. 2.- 



(c) Same as Fig. 2a, but for different values of the total explosion energy E . 




Fig. 2. — (d) Same as Fig. 2a, but for different values of the equipartition factor q. 




Fig. 2. — (e) Same as Fig. 2a, but for different values of the baryon-loading factor To. Note the strong 
dependence of the duration, peak power, and photon energy of the vL v peak on IV 
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Fig. 3. — Dependence of the observables Ho, td and £q on the baryondoading parameter To, the total 
explosion energy Eq = 10 54 /?54 ergs and the external density no (cm~ 3 ), for To ranging from 100 to 1000. 
Other parameters are given in Table 2. 



